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Abstract. The reaction jp — > pyr^r; has been studied with the CBELSA detector at the tagged photon beam 
of the Bonn electron stretcher facility. The reaction shows contributions from zi "'"(1232)77, V(1535)"'"7r'' and 
pao(980) as intermediate states. A partial wave analysis suggests that the reaction proceeds via formation 
of six A resonances, zi(1600)P33, Z\(1920)P33, /\(1700)L»33, ^(1940)Li33, .4(1905)^35, zi(2360)D33, and 
two nucleon resonances A'"(1880)Pii and V(2200)Pi3, for which pole positions and decay branching ratios 
are given. 

PACS: 13. 30. -a Decays of baryons, 13.60.Le Meson production 14.20.Gk Baryon resonances with S = 



1 Introduction 

The study of baryon resonances has found renewed in- 
terest as evidenced by new experimental facilities provid- 
ing continuous photon beams like Jefferson Lab, SPring-S, 
MAMI at Mainz, and ELS A at Bonn. For nucleon res- 
onances up to « 1.7GeV/c^, their abundance and most 
masses are reasonably well explained in quark models even 
though the models use very different assumptions on the 
nature of long-range quark-quark interactions, like effec- 
tive gluon exchange [1] , instanton induced interactions [2] , 
or exchange of Goldstone bosons [3]. At higher masses, the 
models differ in important details of their predictions but 
agree in predicting the existence of many more states than 
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have been found experimentally. Diquark models [4] have 
a reduced number of degrees of freedom and expect fewer 
states [5,6]. Models based on a conformal approximation 
of QCD exploit the correspondence [7] between string the- 
ories on an Anti-de-Sitter (AdS) space and gauge theories 
on its space-time boundary [8,9]. Apparently, AdS/QCD 
predicts a smaller number of states [10-13] but, to our 
knowledge, the problem of missing resonances, of reso- 
nances predicted in quark models but not in AdS/QCD, 
has not yet been discussed. Lattice gauge calculations 
aim at simulating full QCD; substantial progress has been 
achieved [14], however, the calculations do not yet help to 
validate or reject specific quark models. 

In the past, knowledge on nucleon resonances was de- 
rived mostly from partial wave analyses of elastic TriV 
scattering experiments. This method exploits the coupling 
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Fig. 1. Experimental setup at ELS A in Bonn: Electrons delivered by the accelerator ELS A enter the experimental area from 
the left. They hit a radiator and produce a wide-band photon beam. The photon energy is tagged by measuring the electron 
deflecting angle in a dipole magnet. The photons hit protons of a liquid H2 target. Particles, in particular photons, emerging 
from the reaction are detected in the Crystal barrel detector. A scintillation fiber (scifi) detector measures the direction of 
charged particles. A total absorption lead mineral oil Cerenkov counter installed further downstream to monitor the photon 
beam intensity is not shown in the figure. 



of resonances to the Ntt channel; resonances for which 
this coupling is weak are not detectable. The Particle 
Data Listings [15] teach us that at high pion momenta, 
above 2 GeV, high-spin states are formed while formation 
of high-mass states with low angular momenta seems to be 
suppressed. Due to this limitation, only few resonances are 
known so far. Even worse, the evidence for a large fraction 
of them has recently been questioned in a careful analysis 
[16] of a large body of Ntt elastic and charge-exchange 
scattering data. 

The new facilities give access to photo-induced pro- 
duction of resonances; detectors with large solid angle al- 
low one to study decays of resonances into complex final 
states. Thus Ntt is avoided in the entrance and in the 
exit channel. Quark model calculations [17] and first re- 
sults [18,19] seem to support our conjecture that in pho- 
toproduction, excitation of radially excited states might 
be preferred over excitations of high-angular momentum 
states. 

In this paper we present a study of the reaction 

7P pn'^T] (1) 

for photon energies covering the resonance region from 
the pTT^rj production threshold a,t W = 1.63 GeV/c^ up to 
2.5 GeV/c^. In this reaction, decays of A resonances into 
Z\(1232)77 can be studied: the rj acts as an isospin filter 
and resonances decaying into Z\(1232)77 must have isospin 
/ = 3/2. For low photon energies phase space is limited, 
and Z\(1232) and rj should be in a relative S-wave. We thus 
may expect a high sensitivity for baryon resonances with 
isospin 7 = 3/2 and spin J = 3/2 and negative parity. 
If such resonances decay into Ntt, they need L = 2 be- 
tween TV and tt; resonances with these quantum numbers 



are characterized by ^27,2,/ = The lowest mass res- 
onance with these quantum numbers is Z\(1700)L'33. The 
possibility that it may couple to the Z\(1232)77 channel 
was already discussed by Nefkens [20]; his conjecture was 
confirmed recently by the GRAAL collaboration [21]. The 
aim of this paper is to report evidence for A resonances 
above zi(1700)£'33, in particular for a = 3/2^ parity 
doublet consisting of the two resonances Z\(1920)P33 and 
A{1M0)D33. 

2 Experimental setup and data reconstruction 

The experiment was carried out at the tagged photon 
beam of the ELcctron Stretcher Accelerator ELSA at 
Bonn [22] , using the Crystal Barrel detector [23] . A short 
description of the experiment, data reconstruction and 
analysis methods can be found in a recent letter on 27r'' 
photoproduction [24]; further details are given in [25]. A 
schematic drawing of the experimental setup is shown in 
Fig. 1. For the data presented here, ELSA delivered - via 
slow extraction - a continuous electron beam of 3.2 GeV. 
The electrons hit a radiator target of 0.003 (i.e. radi- 
ation length) thickness. 

The energies of photons produced in the radiator tar- 
get were tagged in energy by a measurement of the de- 
flection angle of scattered electrons passing the field of 
a dipole magnet. For singly scattered electrons, the pho- 
ton energy is Ej = Eq — E^- . The tagging detector con- 
sisted of 2 Multi-Wire Proportional Chambers (MWPCs) 
and 14 plastic scintillator bars. The position resolution of 
the MWPCs determined the photon energy resolution of 
0.5 MeV at the highest E^ and 30MeV at the lowest E^, 
whereas the scintillation counters enabled fast timing. The 
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Fig. 2. Left: The distribution of 77 invariant masses after a kinematic fit to 7p — > p47 (with c.l.> 0.1%). There are six entries 
per event. The distribution shows a large peak due to the reaction 7p ptt^tt'^ and two smaller peaks due to 7p — > pn'^r). 
Right: The 77 invariant mass distribution after a kinematic fit to 7p p7r''77 (with c.l.> 10%) and a cut rejecting 7p ptt^tt" 
events (with c.l.> 1%). In the 460-620 MeV/c^ mass interval there are 18589 events, 1120 of them are below the background 
line derived from a linear fit. The fit to 7p prv^r] (with c.l.> 1%) accepts 18379 events and rejects 210 background events 
leading to 910 background events (grey-shaded area). 



energy calibration of the tagger was performed by direct 
injection of the electron beam. 

The photon beam hit a hquid H2 target (length: I = 
52.84mm, diameter: d = 30mm); charged reaction prod- 
ucts were detected in an inner three-layer scintillating- 
fiber detector surrounding the target and positioned at 
mean radii of 5.81cm, 6.17cm and 6.45cm, respectively. 
The 2mm fibers are partly bent to helical shapes (—25°, 
-1-25°, 0° from the inner to the outer layer) to provide an 
unambiguous impact point when a charged particle crosses 
the detector [26]. 

The Crystal-Barrel calorimeter consisted of 1380 16- 
radiation-length CsI(Tl) crystals. The crystals are of trape- 
zoidal shape and point to the center of the target; they 
provide an excellent photon detection efficiency and a high 
granularity. The setup covered about 98% of An (full (f> 
coverage, 12° < 9 < 168°). Its large solid-angle coverage 
allowed for reconstruction of multi-photon final states. 

The first-level trigger was derived from a coincidence 
between the tagging system and the fiber detector. In 
the second-level trigger, a FAst Cluster Encoder (FACE) 
based on cellular logic, provided the number of charged 
and neutral particles detected in the Crystal Barrel. Data 
were taken triggering on events with two (partly three) 
or more particles in the cluster logic. A segmented total- 
absorption oil Cerenkov counter determined the total pho- 
ton flux. 

Events due to reaction (1) were selected by requir- 
ing five clusters of energy deposits in the Crystal Bar- 
rel calorimeter. One of the clusters was required to match 
with the charged particle emerging from the liquid II2 tar- 
get and hitting the scintillation fiber detector. The latter 
cluster was identified as proton, the other four as pho- 
tons. The proton and the four photons were assumed to 



be produced in the target center. More details on event 
reconstruction can be found in [25]. 

These events were subjected to a kinematic fit to the 
7P pAj hypothesis imposing energy and momentum 
conservation. The distribution of the 77 invariant mass of 
one pair against the second pair of surviving events (with 
6 entries per event) is shown in Fig. 2, left panel. Then, the 
7P — *■ p 7r''77 hypothesis was tested and events with a con- 
fidence level (c.l.) exceeding 10% were retained. In a next 
step, events compatible at a c.l. >1% with the 7p — > p2Tr^ 
hypothesis were rejected. The resulting 77 invariant mass 
of the second photon-pair is shown in Fig. 2, right panel. 
These events passed a final kinematic fit to the 7p — > p 7r'^?7 
hypothesis requiring a probability exceeding 1%. The fi- 
nal event sample contains 17469 events due to reaction 
(1) and 910 background events. Events which are likely 
due to the background are identified (and subtracted) by 
selecting those 910 events which are closest in phase space 
to the events falling into the 77 side bins (380-440 MeV/c^; 
640-700MeV/c2). 



3 Mass and angular distributions 

Fig. 3 shows a series of pTr^yy Dalitz plots for M^(prj) and 
M^(7r°?7) versus M^(p7r°) restricted to slices in the "fp in- 
variant mass W ^ with W covering the threshold region (a) 
1.7 <W < 1.9GeV/c2, (b) for 1.9 < I^ < 2.05GcV/c2, 
(c) 2.05 < W < 2.25GeV/c2, and (d) 2.25 < W < 
2.5 GeV/c^. Figs, (e) and (f) show Dalitz plots for A'P{Tn]) 
versus M^(p7r°) for the highest energy bin. The figures in 
the lower part, (g) to (1), show the same Dalitz plots as 
above but calculated from the partial wave fit. All sig- 
nificant structures in the Dalitz plots are reasonably well 
reproduced. At low energies, for W < 2.05GeV/c^, the 
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Fig. 3. Dalitz plot for the reaction piv'^ri for various ranges of the total energy, data (a-f). The plots (a-d) show M^{pri) 

versus M^{pTr''), (e,f) M^{-K^r}) versus lvP{p-K^). With increasing energy, A{12'i2)rj and A'^(1535)7r production fill different 
kinematical regions and are well separated. Af(1535)7r is visible only for high photon energies even though the A'^(1535)7r 
production threshold (~ 1.0 GeV) is lower than the A(12?i2)r] production threshold (~ 1.2 GeV). The lower part (g-f) show 
results of the event-based likelihood fit. 
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Fig. 4. Mass distributions for three values of W = M(7p). Data are represented by dots. The distributions are not corrected 
for the detector acceptance; errors are of statistical nature only. The solid curves show the result of the best fit. The dashed 
line stands for the A{1232)rj, the dashed-dotted line for the A''(1535)5'ii7r, and the dotted line for the pao(980) partial wave 
contribution as determined in the partial wave analysis. 




COS 9(71°) COS 9(r|) cos e(p) 




cos 9(pr|) cos Q{pn°) cos Q{n°r\) 



Fig. 5. Mass and angular distributions for the full data set. The distributions are not corrected for acceptance to allow a fair 
comparison of the fit with the data without introducing any model dependence by extrapolating e.g. over acceptance holes. 
Data are represented by dots, the fit as solid line. Errors are of statistical nature only. The dashed line stands for the D33, the 
dashed-dotted line for the P33 partial wave contribution. There are further PWA contributions which are not shown, (a) prj, 
(b) pir", (c) TT^r] mass distributions. In d)-i) cos 6 distributions are shown. In (d-f), 6 is the angle of a vr" (d), 77 (e), p (f) with 
respect to the incoming photon in the center-of-mass-system (cms); in (g) the angle between the rj and the pr} rest frame, 
in (h) the angle between tt" and 77 in the pn'^ rest frame, and in (i), the angle between tt" and p in the rest frame. 
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Fig. 6. a: Angular distribution of rj with respect to the photon direction in the jp center-of-mass system, b: Center-of-mass 
angular distribution of ao(980) recoiling against the proton with respect to the photon beam axis for E~f > 2.0 GeV. c and d: 
Angular distribution of pions recoiling in the helicity (c) and Godfrey Jackson (d) frame (after side bin subtraction). The latter 
two distributions are presented for 1.9 < < 2 GeV and 1.15 < Mp^ < 1.30 to select Arj decays, the cuts 1.075 < Mp^^ < 1.150 
and 1.300 < Afp^ < 1.375 GeV/c^ define the side bins. Losses due to the cuts are corrected using Monte Carlo simulations (see 
text). The energy region contains large contributions from Z\(1920)P33 and Z\(1940)D33. In all four sub-figures, dots with error 
bars give data, the histogram represents the PWA fit. 



Dalitz plots arc dominated by Z\(1232)r/; at high ener- 
gies, the A^(1535)7r intermediate state is observed. Under 
these kinematical conditions, the two contributions pop- 
ulate different kinematical regions and can be separated 
easily. The early onset of Z\(1232)?7 is, at the first glance, 
surprising since the threshold for Z\(1232)7; is higher than 
for A^(1535)7r production. Within the partial wave analysis 
this is explained by contributions from Z\(1700)I?33 and 
Zi(1600)P33. The two resonances can be produced by an 
electric or magnetic dipole transition (and higher-order 
transitions), respectively, and decay in Z\(1232)?7 via S- 
or P-wave while decays into A^(1535)7r require P- and D- 
waves. 

Three different 2-body invariant masses can be formed 
in reaction (1), the pr], pir'^, and tt^t] mass distributions. 
They are shown in Fig. 4 for three ranges of photon ener- 
gies. The distributions undergo significant changes when 
the photon energy is varied. In the threshold region, (1.7 < 
W < 1.9GeV/c^), the p7r° invariant mass distribution 
(Fig. 4b) suggests a significant zi(1232) contribution (a 
presumption confirmed in the partial wave analysis), the 
prj (Fig. 4a) and n'^ri (Fig. 4c) mass distributions show 
no significant structures. In the 1.9 to 2.1GcV/c^ mass 
region, Z\(1232)7y is still significant, a small prj threshold 
eSect is observed (Fig. 4d). Above 2.1GeV/c2, iV(1535)7r 
- with A^(1535) decaying into Nrj - becomes visible (Fig. 
4g). At these energies, A^(1535) production is significant 
and comparable in intensity to the Z\(1232). In the Tr^yy 
mass distribution, a peak is observed at 1 GeV/c^ (Fig. 



4i) which we assign to the pao(980) intermediate state. 
Mass projections for the full data set and some angular 
distributions are presented in Fig. 5. The P33 and D33 
partial wave contribution are also shown. 

We now turn to the discussion of a few specific angular 
distributions. In the threshold region for pTr'^r] photopro- 
duction we may expect 5-and P-waves to dominate the 
decay mode. S'-wavc decays into Z\(1232)77 require 
quantum numbers which can be reached by an electric 
dipole transitions; magnetic dipole transitions lead to P33 
quantum numbers for which P-wave decays into Zi(1232)77 
are possible. If one of the two amplitudes prevails, we may 
expect an isotropic angular distribution for the recoiling 
pion. In case of both amplitudes being present, S-P inter- 
ference leads to a linear increase with cos 6rj , a distribution 
which is clearly contributing to the data (Fig. 6a). 

At the largest photon energies, a peak due to ao(980) 
production can be identified in the irrj invariant mass dis- 
tribution. In Fig. 6b the ao(980) angular distribution with 
respect to the photon beam is presented (after side bin 
subtraction). The ao(980) signal region is defined by 0.96 < 
A'U,, < 1.04GeV/c2, the sidebins with 0.88 < M„j < 0.96 
and 1.04 < Af^j, < 1.12 are subtracted with weight 1/2. 
The side bin subtraction suffers from the large ao(980) 
width, hence the distribution must be viewed with some 
precaution. However, it agrees reasonably well with the 
PWA result and the subtraction is justified a posteriori. 

Surprisingly, ao(980) is not mostly produced in forward 
direction as might be expected if the main ao(980) produc- 
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tion mechanism would have been vector-meson exchange 
in the i-channel (exploiting a ao(980) j p/uj coupling). 
A significant fraction of the ao(980) signal seems to stem 
from the decay of high-mass resonances. This suspicion is 
supported by the partial wave analysis. 

Fig. 6c and d show Z\(1232) — > p7r° decay angular dis- 
tributions which should reflect the spin J = 3/2+ quan- 
tum numbers of zi(1232). The tt'^ angular distribution 
from A decays is shown as pion emission angle with re- 
spect to the A flight direction (helicity frame, Fig. 6c), 
or with respect to the photon beam direction (Godfrey- 
Jackson frame, Fig. 6d). The photon energy is restricted 
to W < 1.9GeV/c^. The distributions are essentially flat 
(preventing a straightforward identification of the Z\(1232) 
spin) indicating that Z\(1232) are likely not produced by 
one helicity amplitude only and in a single reaction chain. 
Explicit expressions for typical angular distributions arc 
given elsewhere [28]. 

4 Partial wave analysis 
4.1 The method 

The qualitative findings discussed in the previous section 
are confirmed in a partial wave analysis (PWA). The for- 
malism is documented in [27-29]. In addition to the data 
presented here, data on photoproduction of single pions, 
?7's, and of A and S hypcrons are included as well as some 
important partial-wave amplitudes for ttN elastic scatter- 
ing. References to the data, an outline of the PWA method 
and the definition of total likelihood and likelihood con- 
tributions can be found in [30]. The different data sets 
entered with appropriately chosen weights Wi ; the weights 
vary from 1 (for data which should be reproduced qualita- 
tively) to 30 (for low-statistics data like polarization data 
which we insist to be well reproduced in the fits). The 
weights are given in Table 4 of [30] . The weights are chosen 
to ensure that there is good qualitative agreement between 
data and fit even for low-statistics data which otherwise 
may be dominated by the influence of high-statistics data. 
The mean weight is W = 4.2. For this data, w^p_p ^o^ = 10 
was chosen; (much) smaller weights lead to discrepancies 
between fit and this data, higher weights deteriorate the 
fits elsewhere. We minimized a pseudo-log-likelihood func- 
tion In £tot defined as a sum of the log-likelihoods of fits to 
7P — *■ tt^tt" and jp — > 7r°?7, and of the x^/2 contributions 
of all data given in the form of histograms. 

The 7p — > pn'^TT^ and jp — > pir^rj data are fitted using 
an event-based maximum likelihood method which takes 
into account the full correlations between all variables in 
the 5-dimensional phase space. New data on the beam 
asymmetry S for ^p — > pTr^rj [31] were included in the 
partial wave analysis. Since the prediction for S were al- 
ready close to the data [32], the new data on E did not 
change the PWA results. 

As a starting point, we used the set of resonances of 
[30] and added resonances in the partial waves listed in 
Table 1. The Table gives the quantum numbers of baryon 
resonances which could contribute to the pirr] final state 



Table 1. The lowest partial waves which may contribute to 
7P ^ p ivq. 



L A{1232)r] Af(1535)7r pao(980) 

D33 Pii,Pii Pii, P31 

1 P31, P33, F35 Sii, D13, S31, -D33 Sii, D13, S31, D33 

2 S31, D33, -035, G37 Pl3, Fi5, P33, F35 Pl3, Fis, P33, F35 



assuming that the relative orbital angular momenta be- 
tween meson and baryon in the three visible intermediate 
states Z\(1232)?7, A^(1535)7r, and pao(980) are restricted 
to i < 2. First, resonances of the first line were tested 
but their decays with L = 1 and L = 2 were also ad- 
mitted. It turned out that for nearly all resonances, in- 
clusion of L = 2 decays led to marginal likelihood im- 
provements (ziln/^TrOj) < 5 corresponding to a statisti- 
cal significance of less than two standard deviations. The 
Z\(1232)?7 D-wave contribution, e.g., was found to be negli- 
gible for A{1700)D33, below 1% for Z\(1940)i:»33 while for 
Z\(2360)D33, the S and D wave couplings had similar mag- 
nitudes. The partial waves in the second (and third) line of 
Table 1 were tested one by one. No significant contribution 
was found except for a small fraction of Z\(1905)F35 in its 
Z\(1232)7; decay. Its fraction in iV(1535)7r (with L = 2) is 
marginal. 

The quality of the final PWA fit to the data is shown 
in Fig. 4 and 5. In the latter figure, the full statistics inte- 
grated over the full energy range is shown; the agreement 
between data and fit result is very convincing and equally 
good when mass slices are selected or other variables are 
plotted. It has to be remembered that a likelihood fit 
uses the full kinematics of each event and that neither 
the Dalitz plots (Fig. 3) nor the projections (Fig. 4) are 
fitted. 



4.2 Isospin considerations 

Some resonances, even though known to exist, cannot be 
reasonably included in the partial wave analysis, for tech- 
nical reasons. We need e.g. a = 1/2+ resonance at 
about 1900MeV/c^ A iV(1880)Pii resonance (formerly 
called iV(1840) by us [19]) is required to describe the two 
reactions ^p ^ A and 7p — > 17 (which define the 
isospin of contributing resonances). Both, A^(1880)Pii and 
Z\(1910)P3i, can decay into iV(1535)7r and/or pao(980), 
and there is no information in these two decay modes to 
decide on the isospin. Hence both might be present. The 
Z\(1910)P3i might decay into Z\(1232)?7 but this is not 
observed. Thus there is evidence for A^(1880)Pii and no 
direct evidence for an additional Zi(1910)P3i. However, 
a Z\(1910)P3i might nevertheless exist and contribute to 
p7r°?7 if it couples to A^(1535)7r and/or pao(980) but not 
to S and Z\(1232)?7 (or only with a weak coupling). 
In such a case (which admittedly is a bit constructed but 
not forbidden by any selection rule), Z\(1910)P3i decays 
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Fig. 7. Total cross sections for pix rj. In both panels, our data are represented by • with error bars; the solid curves 

represents our PWA fit. Comparison to other data and two decompositions of the total cross sections are shown in two sub- 
figures, a) Our data compared to data from Nakabayashi et al. [33] ■. Curve 1 shows the cross section for yp —> A{1232) r], 
curve 2 that for A'^(1535)7r, and curve 3 that for 7p — > pao(980). These cross sections are not corrected for unseen decay modes, 
b) open circles: Ajaka et al. [21]. Partial wave decomposition of the cross section into the D^s (curve 4) and the P33 partial 
wave (curve 5). 



into pir^i] are discarded. The reason for this is that in- 
troduction of two close-by resonances in one final state 
having different isospin but otherwise identical quantum 
numbers leads very often to very large amplitudes for both 
resonances and destructive interference. 

The 7V(2200)Pi3 was introduced, too, in [19]. Here, it 
is found to contribute significantly to pir^rj but not to 
Z\(1232)7^. Hence we assume it to belong to the familiy of 
nuclcon resonances, and do no include a further Z\(2200)P33 
For the same reason we do not allow A^(1875)£'i3 to de- 
cay into A^(1535)7r and/or pao(980); in the 3/2~ partial 
wave, the intensity inp7r'^?7 is assigned to Z\(1940)I?33. The 
isospin / = 3/2 is identified from the Arj decay mode. Con- 
tributions from / = 1/2 to A^(1535)7r and/or pag (980) are 
possible but are not considered here because of limitations 
of the data base. A separation of two isospin contributions 
would require a second isospin channel, e.g., the study of 
7P — > nir^r]. 



4.3 The total cross section 

The final fit required contributions to the p — > pir^r] final 
state from eight resonances; these are listed in Table 2. 
Two of them do not couple to zi(1232)?7 and are inter- 



Table 2. Resonances used in the partial wave analysis. 



iV(1880)Pii 


Ar(2200)Pi3 


Z\(1600)P33 


Z\(1920)P33 


Zi(1700)7?33 


Z\(1940)D33 


Zi(2360)D33 


Z\(1905)P35 



preted as nuclcon resonances. The Af(1880) could be the 
missing Pu resonance forming, jointly with A^(1900)Pi3, 
iV(2000)Fi5, iV(1990)Fi7, a super- mult iplet with (domi- 
nantly) L=2, 5=3/2. The Af(2200)Fi3 might be a iV(1900)Pi3 
radial excitation. 

In Fig. 7 the total cross section is displayed. The points 
with errors give the acceptance-corrected results of the 
measurement and their statistical and systematic errors 
but not the ±15% systematic error assigned to the photon- 
fiux normalization. The solid curve shows the result of 
the partial wave analysis (PWA). The total cross section 
reaches the maximum of 4 /Ltb in the 2 GeV/c^ region and 
then decreases slowly to 3 fih. In the threshold region, our 
cross section is fully compatible with results obtained at 
Sendai [33]. The cross section determined by GRAAL [21] 
falls systematically below our values but the difference is 
covered by the 15% normalization error. 

The acceptance is calculated using Monte Carlo tech- 
niques. The full apparatus is simulated using GEANT3. 
Detection and reconstruction efficiency, probabilities for 
split-offs for photons and charged particles, confidence 
level distributions for the different kinematical fits for real 
data and Monte Carlo data were compared carefully. The 
detector is well understood by the Monte Carlo simulation 
[25] . The mean acceptance for the pA-f final state is 15-20% 
but the acceptance vanishes for forward protons escaping 
through the forward hole in the detector. Hence a proper 
extrapolation of the four-dimensional event distribution 
into the uncovered region of the detector is essential. 

For this purpose, phase space distributed events due 
reaction (1) were simulated according to the photon flux. 
These events were weighted with the PWA solution. The 
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resulting event sample represents the 'true' physics, in the 
limit that data and partial wave analysis are correct. The 
curves in Figures 6 and 7 are calculated as Monte Carlo 
integral over this (weighted) event sample. The detection 
efficiency for the data is obtained from the ratio of re- 
constructed and generated events of the weighted Monte 
Carlo event sample. 

The systematic error comprises the error in the recon- 
struction efficiency (5.7%) and a few small errors (« 1%) 
due to uncertainties in target position and profile of the 
photon beam, and the extrapolation of the intensity dis- 
tribution into blind regions of the detector (mainly for 
forward protons). The latter uncertainty is estimated to 
about 5% by using different PWA solutions giving accept- 
able descriptions of the data. These errors are added in 
quadrature and included in Fig. 7. The largest error stems 
from the uncertainty in the photon flux: the normaliza- 
tion derived from a comparison of the photoproduction 
cross section with SAID (see [25]). The normalization er- 
ror of ±15% is not included. 

The excitation functions for the three intermediate 
states A{1232)t], 7V(1535)7r, and pao(980), deduced from 
the partial wave analysis, are shown in Fig. 7a. The Z\(1232)?7 
makes the most significant contribution to the total cross 
section; at threshold, it dominates the reaction. It reaches 
a peak height just below 2GeV/c'^. The 7V(1535)7r in- 
termediate state exhibits a structure at 1.9GeV/c^ (due 
to the iV(1880)Fii) and a second bump at 2.2GeV/c2 
which we interpret as 7V(2200)Pi3. With increasing mass, 
pao(980) gains a notable intensity. The excitation func- 
tions shown in Fig. 7 represent the best solution. It should 
be stressed that these contributions arc qualitative; the 
numerical values differ significantly for different fits, much 
more than the pole positions of the contribution reso- 
nances. The errors in the excitation functions can be, in 
particular at high energies, as large as 30%. 

4.4 Resonances contributing to p-7v°r] 

In a first approach, we used Breit-Wigner resonances for 
less important waves while Su, Pn, P13, i-33, and D33 
were described within the K-matrix/P-vector approach 
[30]. For these waves, elastic ttN scattering amplitudes 
from [16] were included in these fits. For the two waves 
P33 and £'33 discussed here, a satisfactory description was 
obtained with a 6-channel (A''7r, Z\(1232)7r (5, -D-waves 
for D33 and P, P- waves for P33), Z\(1232)?7, iV(1535)7r, 
A'^ao(980)) 3-pole parameterization. Parameterizations of 
the Sii, Pii, and P13 amplitudes are given in [29], those 
for P33 and are presented in Tables 3 and 4. 

All eight resonances are required to get a good fit. Ta- 
ble 5 lists the resonances used in the final fits (column 1 
and 2), their pole positions (column 3) and their decay 
modes leading to the pn^r] final states (column 4). The 
fractional contributions (adding up to 100% of the inte- 
grated p7r°?7 cross section) are listed in column 5, while 
columns 6 and 7 give the likelihood change, for the to- 
tal data set and for the data presented here, respectively, 
when a resonance is assumed not to couple to pir^rj. 



Table 3. Properties of the P33 Tv-matrix poles. The masses 
Mi and if-matrix coupling constants Qa^ from resonance (i) to 
the decay mode a are given in GeV, the helicities A in GeV~ 2 , 
s in GeV^/c*; the / - describing N-k — > A-k contact interactions 
- are dimensionless. The parameters marked with * were fixed 
in the fit. 







Pole 1 


Pole 2 


Pole 3 






h'h 


1232±3 


1795±40 


2020±80 






Ax/2 


-0.125 


0.031 


0.026 








-0.264 


0.123 


0.081 




a 




ya 


(2) 
Qa 


(3) 

ga 


/la 


1 


iV(940)7r 


1.20 


1.25 


0.92 


-1.40 


2 


zi(1232)7r(P) 


-0.52 


-1.73 


-0.69 


0.81 


3 


Z\(1232)7r(F) 


0.* 


0.19 


-0.76 


0.* 


4 


zi(1232)7) 


0.* 


0.75 


-0.63 


0.* 


5 


iV(1535)7r 


0.* 


0.94 


1.04 


0.* 


6 


iV(940)p 


0.* 


0.05 


-0.06 


0.* 


Table 4. Properties of the Dyj, /sT-matrix poles. 






Pole 1 


Pole 2 


Pole 3 






Mi 


1780±50 


1980±60 


2400±90 








-H0.041 


0.205 


0.228 






A3/2 


0.252 


-0.003 


0.143 




a 






ya 


ya 


/la 


1 


7V(940)7r 


0.31 


0.66 


1.22 


-0.62 


2 


Zl(1232)7r(5') 


0.51 


-0.72 


-0.16 


-0.17 


3 


Zi(1232)7r(I>) 


-1.21 


-1.49 


-1.51 


-0.63 


4 


Zl(1232)?7 


-0.39 


-0.23 


0.81 


0.* 


5 


iV(1535)7r 


-0.38 


-0.40 


-0.78 


0.* 


6 


A/'(940)p 


-0.14 


-0.26 


-0.24 


0.* 



In the P33 wave we find a first pole above Z\(1232) 
at 1510+5° - ill5 ± 20MeV/c2 and a further pole at 
1980;^^ - a75l^^MeV/c2. If the couplings of the first 
pole to p-KT] are reduced to zero, the likelihood for the 
7P pirr] data deteriorates by 60. Removing pirri cou- 
plings of the highest K-matrix pole "costs" 63 units in 
likelihood. If the latter pole is completely removed from 
the fit, the likelihood value for ttt] channel changes by 177 
and total likelihood by 2230. Such fits are unacceptable. 

Using only Zi(1700)P'33 for the wave, no satisfac- 
tory description of the data was obtained; a second res- 
onance proved to be essential. The pole position for the 
second resonance was found at (1985 ± 35) — i (195 ± 25) 
MeV/c^. The likelihood value for the jp pn'^rj reac- 
tion changed by AlnC^p^py^o^ =151, the total likelihood 
by 2005 units; these are highly significant numbers. With 
Z\(1940)i:>33 and without Z\(1700) D33, the fit was also un- 
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Fig. 8. The nN elastic scattering amplitude in the P33 (left) and -D33 (right) wave [16] and our fit. The P33 shows the distinct 
phase motion for Zi(1232); a second resonance at 1600MeV/c^ is difficult to see, and there is no sign for a third resonance. 
Likewise, there is a rapid phase motion at 1700MeV/c^ in the D33 but no sign for any resonance above it. 

Table 5. Resonances contributing to yp — > pirrj, their pole positions and decay modes, and change of likelihood when resonance 
couplings to Nivr] are reduced to zero. A change in likelihood by 50 corresponds to a change of 25 or to a statistical significance 
of 5 a. 



Name 


Wave 


Pole (MeV/c^) 


Decay modes 


Fraction 


A In Ctot 


Z\ln£^ 


7V(1880) 


Pii 


1880 - illO 


pao,N{1535)Tv 


12% 


2050 


70 


^■(2200) 


Pl3 


2200 - il50 


pao,iV(1535)7r 


7% 


1392 


73 


Z\(1600) 


P33 


1510 - ill5 


Ari 


1 22% 


884 


60 


Z\(1920) 


P33 


1980 - il75 


Ar],pao 


1818 


63 


Z\(1700) 


1533 


1640 - il60 


Ari, N{1535)tv 




2611 


84 


Z\(1940) 


1533 


1985 - il95 


Ar],pao,N{1535) 


1 24% 


1124 


52 


Zi(2360) 


1533 


2320 - i250 


Arj,pao,N{1535) 




1018 


24 


zi(1905) 


F35 


1920 - il45 


Ari,N{1535)-K 


7% 


1991 


65 




f-channel 


p-io-, ti-channel p- 


Z\- exchange 


28% 


6340 


208 



acceptable (ziln £^p_>p^o^ =137, ZilnAot =4782). The 
fit improved notably (Z\ln C^p^pT^Oj^ =28, Z\ln£tot=770) 
by adding a third D33 state. Its mass optimized at about 
2.35GeV/c2 and a width in the 400-600 MeV/c^ range. Its 
mass is close to the end of the available phase space. An 
upper bound for its mass is not well defined; hence we do 
not claim that this state exists. 

When an additional iV(2000)Di3 resonance decaying 
into iV(1535)7r and A^ao(980) was introduced, strong in- 
terference between the two amplitudes was observed which 
increases the systematical errors of the PWA results. A 
precise mass determination of the Z\(1940)D33 will hence 
require separation of the two isotopic resonances, e.g. by 
measurements of •jN n-K^-q. This ambiguity affects of 
course the branching ratios as well. 

We assign the absence of Z\(1920)P33 and Z\(1940)D33 
in the phase shift analysis of Arndt et al. [16] to their weak 
couplings to the Ntt channel. To check if our data are 
consistent with elastic scattering amplitudes, we included 
these amplitude in the fit. The result is shown in Fig. 8: 
the existence of the two resonances is not in conflict with 
elastic scattering. 

Fig. 3e shows a structure in the pir^ invariant mass 
at about 1.5 GeV (while Fig. 3f would suggest a higher 



mass) . We tested the possibility that the structure signals 
contributions from jp — A^(1520)77 or jp — > Zi(1600)77. 
Both attempts yielded marginal improvements in likeli- 
hood and contributions of about 0.5%. The PWA inter- 
prets the structure as interference between pao(980) and 
7V(1535)7r. 

To check the stability of the solution we added, one 
by one, new resonances to the analysis. Adding a reso- 
nance with 5*31 or 6*37 quantum numbers returned masses 
at the end of the phase space, very broad widths, negli- 
gible fractional contributions and marginal improvements 
in likelihood. For a D35 wave, stable fit results were ob- 
tained with M w 2150 MeV, T « 160 McV, and a 1% 
fractional contribution. The improvement of the likelihood 
was marginal, < 10. Adding a fourth pole to the D33 wave, 
reproduced the T-matrix pole position reported above and 
created an additional pole in the region of 2 GeV/c^. How- 
ever the statistical evidence for this was again marginal, 
likelihood changed by < 10. The present data do not sup- 
port the need to introduce additional resonances. We tried 
to improve the fits by adding decays of Z\(2360)D33 into 
Z\(1600)P33^ or iV(2200)Pi3 iV(1710)?7. The statistical 
significance for these transition was weak, and the frac- 
tional contribution to the data stayed below 0.5%. 
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M(P33) [MeV/c^] M(P33) [MeV/c^] 

Fig. 9. Likelihood as a function of the mass of a resonance. The width is fixed to 350MeV/c^. The left plot gives the 

change of the total likelihood, the right plot the change of the likelihood contribution from jp pvr"r;, from this data and from 
beam asymmetry data [31]. The minimum (reference) value is at about -280 000 (a) and -6 000, respectively. 




M(D33) [MeV/c^] M(D33) [MeV/c^] 



Fig. 10. Likelihood as a function of the mass of a AD33 resonance. The width is fixed to 350MeV/c^. The left plot gives the 
change of the total likelihood, the right plot the change of the likelihood contribution from 7p — > pi""'), this data and beam 
asymmetry data [31]. 



4.5 Scans 

In a second analysis we replaced, alternatively, the P33 or 
D33 partial wave by a sum of a K-matrix for the energy 
region below 1.85GeV/c^ and a relativistic Breit-Wigner 
(BW) amplitude describing the higher-mass region. The 
mass of the P33 Breit-Wigner state was found to be 1995 ± 
50 MeV and the width 380 ± 60 MeV. The BW mass of 
the second D33 resonance was found to be 1970 ±50 MeV 
and the width 350 ± 70 MeV. 

The findings can be visualized in scans showing the 
changes of — In Ctot obtained in the fit where the mass of 
a resonance under study is changed in steps and fixed 
while all other variables are allowed to vary freely. In 
Fig. 9, the P33 amplitude is described by three BW reso- 
nances, Z\(1232). Z\(1600) and a third resonance the mass 
of which is scanned. The mass of this third resonance was 
scanned from 1.85 to 2.5GeV/c^ and its width fixed to 
350MeV/c^. Fig. 9 (left) shows the change in total hke- 
lihood. Fig. 9 (right) the likelihood chance for the data 
on 7p — > pir^r]. A clear minimum - which we identify 
with Z\(1920)P33 - is observed in both cases; the data on 
7P pn^rj have a preference for a somewhat higher mass. 

In Fig. 10, a second D33 resonance is added to zi(1700). 
Again, the mass scan covered the 1.85 to 2.5GeV/c^ re- 
gion, and the width was fixed to 350MeV/c^. Wide min- 
ima are observed in the 1900-2020 MeV/c^ region. The 



1)33 wave shows a second minimum at about 2.4GeV/c^, 
in agreement with the fits described above. The first min- 
imum is interpreted as evidence for Zi(1940)£'33, the sec- 
ond one may indicate contributions from a higher mass 
resonance in this partial wave. 



4.6 PWA results 

Table 6 summarizes our results. The Breit-Wigner masses 
and widths and the pole positions are mostly compatible 
with Review of Particle Properties (RPP) [15]. Most de- 
cay branching ratios are reported here for the first time. 
They are calculated as ratio of partial decay width and 
total width at the position of the Breit-Wigner mass. The 
best known states in Table 6 are z4(1600)P33, Z\(1700)I?33, 
and Z\(1920)P33, all 3-star or 4-star resonances in [15]. The 
high rating does however not correspond to a good agree- 
ment in resonance properties derived in different analy- 
ses. In Table 7 we compare our pole positions with those 
of Hohler et al. [34], Cutkosky et al. [35], and Arndt et 
al. [16]. It is of course extremely intriguing that three out 
of five of them were not observed in the most recent anal- 
ysis [16]. The Z\(1600)P33 iiT-matrix coupling to Ar] and 
iV(1535)7r are sizable (see Table 3); when the pole is ap- 
proximated by a Breit-Wigner resonance, its couplings to 
these two channels become very small due to phase space 
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Table 6. Pole positions, Breit-Wigner parameters, and decay branching ratios of A resonances. Masses and widths are given 
in MeV/c^, branching ratios in %. The branching ratios are corrected for Clebsch-Gordan coefficients and unseen decay modes 
of final-state mesons (tt" and rj). The helicity couplings are in GeV"^''^. 





^'ipole. 


rpole 


Mbw 


Y^BW 
^ tot 


Brjv,r 




Brjv(1535)-,r 


Brivc 


io(980) 


Ai/2 


A3/2 


zi(1600)P33 


15101^° 


230 ± 40 


1650 ± 40 


530 ± 60 


10 ±3 


0" 












not well defined 


zi(1920)P33 


1980145 


sioi^S 


1990 ± 35 


330 ± 60 


15 ±8 


10 ±5 


6±4 


4 


± 


2 


22 ±8 


42 ± 12 


Zi(1700)D33 


1640 ± 25 


325 ± 35 


1790 ± 30 


580 ± 60 


20 ±7 


2± 1 


4±2 









160 ± 40 


150 ± 30 


Z\(1940)D33 


1985 ± 30 


390 ± 50 


1990 ± 40 


410 ± 70 


9±4 


4±2 


2± 1 


2 


± 


1 


160 ± 40 


110 ±30 


Z\(2360)D33 


2320 ± 60 


500 ± 100 


2310 ± 60 


640 ± 120 


10 ±4 


60 ±20 


7±4 


9 


± 


4 


20 


fa 40 



" The partial width at the Breit-Wigner mass vanishes due to phase space arguments. The integrated contribution is in the 
order of a few %. 



Table 7. Pole positions (in MeV/c^) of "well established" A 
resonances and the one-star Z\(1940)-D33. 



Resonance 


Hohler 


Cutkosky 


Arndt 


this work 


Zi(1600)P33 


1550 


1550 


-ilOO 


1457 - i200 


1510 


- ill5 


Zi(1700)D33 


1651 - i8C 


) 1675 


- illO 


1632 - il26 


1650 


- il60 


zi(1920)P33 


1900 


1900 


-il50 


not seen 


1980 


- il55 


/i(1905)F35 


1829 -il52 1830 


-il40 


not seen 


1920 


- il45 


Z\(1940)D33 


not seen 


1900 


-ilOO 


not seen 


1985 


-il95 



onances is dissipated, into high-momenta of the daughter 
particles, into high-mass mesons or into high-mass daugh- 
ter baryon resonances. The results suggest that high-mass 
parent baryons seem to take advantage of all allowed de- 
cay modes, including intermediate meson and baryon reso- 
nances even at rather high masses. The situation resembles 
NN annihilation [38,39]: in this process, the production 
of high-mass mesons is preferred in comparison to pro- 
duction of low-mass mesons even though the phase space 
is of course larger for the latter reaction. Likewise, there 
seems to be no preference for the kinematically favored 
Nt: decay mode. 



limitations. The K-matrix pole and the Breit-Wigner mass 
differ considerably; this fact renders a definition of the he- 
licity couplings difficult and we refrain from giving num- 
bers. 

The Z\(1920)P33 is the third resonance in this partial 
wave, after the well-known Z\(1232) and the Ropcr-likc 
Zi(1600)P33. Z\(1920)P33 was observed in several analy- 
ses including those of Hohler et al. [34] and Cutkosky 
et al. [35] in elastic ttN scattering data, and by Manley 
et al. [37] in a combined analysis of ttN elastic scatter- 
ing and pion induced double-pion production who found 
M=2057±110, r=460±320MeV/c2. In RPP [15] it is ranked 
as 3-star resonance. In the analysis by Arndt et al. [16], 
Z\(1920)P33 was not found. Z\(1600)P33 is observed in aU 
four analyses of Table 7; in [19], it was not needed to get 
a good description of the data taken into account at that 
time. 

The Z\(1940)L'33 state has only one star in the RPP 
classification. It was observed by Cutkosky et al. [35] in 
elastic ttN scattering data with mass M=1940±100MeV/c^ 
and width P = 200± lOOMeV/c^, and by Chew et al. [36] 
at M=2058 ± 35MeV/c2 and P:=198±46MeV/c2 while 
Manley found M=2071±100MeV/c2,P=460±320McV/c2 
[37]. In [16], Z\(1940)i:>33 was not found. 

The confirmation of the two resonances Z\(1920)P33 
and Z\(1940)P'33 in a new reaction is a major result of this 
analysis. In addition, several partial decay branching ra- 
tios have been determined. These provide first insight into 
an old question on how the mass of high-lying parent res- 



5 Summary and conclusions 



We have reported a study of the photoproduction process 
7P ^ pir'^r] at the Electron Stretcher Apparatus ELSA 
at Bonn. About 16000 pir^r] events were extracted above 
a small background. For small photon energies, the cross 
section for the reaction agrees with previous results from 
Sendai and Graal. In this energy region, the process is 
dominated by the Z\(1232)?7. At higher photon energies, a 
significant fraction of the process proceeds via the iV(1535)r7 
and some pao(980) contribution is observed. 

A partial wave analysis determines the resonant contri- 
butions. Two waves with P33 and D33 quantum numbers, 
with = 3/2+ and 3/2~, dominate the reaction. We find 
evidence for two resonances Z\(1920)P33 and Zl(1940)D33 
and indications for a further high-mass £^33 resonance. 
The evidence for the two states was already communi- 
cated in a letter publication [40] . 

A detailed discussion of the role of zi(1920)P33 and 
Z\(1940)P'33 within the spectrum of A resonances will be 
given elsewhere [41]. Here we just mention the main re- 
sults. 

1. In [41], the 3 negative parity resonances Z\(1900)S'3i, 
Z\(1940)P'33, and Z\(1930)P)35 are interpreted as triplet 
of resonances with L — 1, S — 3/2, N = 1 coupling to 
J = 1/2, 3/2, 5/2 and belonging to the (P>,L^)5 = 
(56, 13^)3/2 supcrmultiplct. Alternatively, Zi(1900)53i 
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and A{1M0)D33 could be a L = 1,5" = 1/2, N = 
1 doublet, and zi(1930)i:'35 and Zi(2200)G'37 a L = 
3,5" = 1/2 doublet. The mass of the triplet is rather 
low compared to quark model predictions [1,2] while 
the Z\(2200)G37 mass is perfectly compatible with the 
models. 

2. The two resonances Z\(1920)P33 and Z\(1940)D33 form 
a parity doublet [42] . Parity doublets are of topical in- 
terests, we refer the reader to two recent reviews [44, 
45]. In the latter review, there is a written request 
to search for the doublet reported here. The conjec- 
ture that the reason for the mass degeneracy is due to 
a restoration of chiral symmetry in high-mass baryon 
states falls of course outside of experimental observ- 
ability, and the interpretation is controversially dis- 
cussed [46-48]. 

3. All masses in the A spectrum are perfectly compatible 
with models based on AdS/QCD when confinement is 
modeled by a soft wall in the holographic variable [12]. 
The hard wall approximation [11] gives a good general 
survey but fails in important details. 
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